Second sound and the superfluid fraction in a resonantly interacting Fermi gas 
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Superfluidity is a macroscopic quantum phenomenon, which shows up below a critical temperature 
[1, 2] and leads to a peculiar behavior of matter, with frictionless flow, the formation of quantized 
vortices, and the quenching of the moment of inertia being intriguing examples. A remarkable 
explanation for many phenomena exhibited by a superfluid at finite temperature can be given in 
terms of a two-fluid mixture [3, 4] comprised of a normal component that behaves like a usual 
fluid and a superfluid component with zero viscosity and zero entropy. Important examples of 
superfluid systems are liquid helium and neutron stars. More recently, ultracold atomic gases have 
emerged as new superfluid systems with unprecedented possibilities to control interactions and 
external confinement. Here we report the first observation of 'second sound' in an ultracold Fermi 
gas with resonant interactions. Second sound is a striking manifestation of the two-component 
nature of a superfluid and corresponds to an entropy wave, where the superfluid and the non- 
superfluid components oscillate in opposite phase, different from ordinary sound ('first sound'), 
where they oscillate in phase. The speed of second sound depends explicitly on the value of the 
superfluid fraction [5], a quantity sensitive to the spectrum of elementary excitations [6]. Our 
measurements allow us to extract the temperature dependence of the superfluid fraction, which in 
strongly interacting quantum gases has been an inaccessible quantity so far. 



Until now, second sound has only been observed in 
superfluid helium [7], the paradigm of quantum fluids 
characterized by strong interactions, for the description 
of which Landau developed his theory of two- fluid hydro- 
dynamics [4] . In liquid helium, second sound can be gen- 
erated [8] by local time-dependent heating and detected 
by observing the propagation of the resulting tempera- 
ture wave. Second sound is an essentially isobaric oscil- 
lation, in contrast to the adiabatic nature of first sound. 

In experiments with ultracold atomic quantum gases, 
the observation of second sound has been a long-standing, 
but so far elusive goal. In dilute Bose-Einstcin con- 
densed samples, the relative motion of the condensate 
with respect to the thermal component represents a re- 
lated effect, the observation of which has been reported 
in refs. [9, 10]. However, a direct observation of second- 
sound waves has not been achieved. In resonantly inter- 
acting Fermi gases [11, 12], superfluidity and the univer- 
sal thermodynamics [13-17] have been subjects of intense 
research. In these systems also the normal component be- 
haves in a deeply hydrodynamic way, which means that 
Landau's two-fluid theory can be readily applied. This 
suggests a behavior similar to superfluid helium, includ- 
ing the occurrence of second sound. 

Our system is an ultracold, superfluid sample of 
fermionic 6 Li atoms, prepared in a highly elongated har- 
monic trapping potential (Methods) by well-established 
procedures of laser and evaporative cooling [18]. The 
sample consists of N = 3.0 x 10 5 atoms in a balanced mix- 



ture of the two lowest spin states, and is about 500 /im 
long and 20 ^m wide. It is characterized by the Fermi 
temperature Tp Rp ss 0.9 [iK (Methods). A magnetic bias 
field of 834 G is applied, which tunes the interaction be- 
tween the two spin components right on top of an s-wave 
scattering length resonance [19] ('unitarity limit' of in- 
teractions). The cloud's temperature T is determined by 
analyzing axial density profiles [16, 18], using the knowl- 
edge of the equation of state (EOS) from ref. [17]. 

Our method for observing sound propagation builds 
on the classical scheme to detect the propagation of first 
sound in Bose-Einstein condensates [20], which was also 
applied to resonantly interacting Fermi gases [21]. The 
general idea is to prepare the quantum gas in a trap that 
is highly elongated, to create a local perturbation, and to 
detect its one-dimensional propagation as a pulse along 
the long trap axis. In the case of first sound, both the cre- 
ation and the detection of such an excitation are straight- 
forward, while being less obvious for second sound. 

For the local excitation of the cloud, we use a repul- 
sive dipole potential created by a tightly focusscd green 
laser beam (Methods) that perpendicularly intercepts the 
trapped sample in its center, as shown in Fig. la. To ex- 
cite first sound, we suddenly turn on the repulsive beam. 
The local reduction of the trapping potential acts on the 
superfluid and normal component in the same way and 
creates a small hump in the axial density distribution. To 
excite second sound, we keep the green beam's power con- 
stant during the whole experimental sequence with the 
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FIG. 1. Observing the propagation of first and sec- 
ond sound, a, The basic geometry of exciting the optically 
trapped cloud with a weak, power-modulated repulsive laser 
beam (green), which perpendicularly intersects the trapping 
beam (red). The trapped cloud has a superfluid core, sur- 
rounded by a normal region, b and c, Normalized differential 
axial density profiles 5ni (z,t)/ni jmax taken for variable de- 
lay times after the excitation show the propagation of first 
sound (local density increase, bright) and second sound (lo- 
cal decrease, dark). The temperature of the atomic cloud is 
T = 0.135(10)T£. rap . The vertical dashed lines indicate the 
axial region where superfluid is expected to exist according 
to a recent determination of the critical temperature [17]. 



exception of a short power- modulation burst, which con- 
tains 8 sinusoidal oscillations in 4.5 ms (Methods). The 
fast modulation locally drives the system out of equilib- 
rium and the following relaxation increases entropy and 
temperature. The duration of the burst is chosen such 
that the system can establish a local thermal equilibrium 
on a length scale that covers the transverse cloud size, but 
is much shorter than the axial extension of the cloud. In 
all cases, we take care that the excitation remains a small 
perturbation of the whole system, which globally stays in 
a thermal equilibrium state. 

For detection we record the axial density profile n\ (z, t) 
for various time delays t after the excitation pulse, where 
n\{z,t) is the number density integrated over the trans- 
verse degrees of freedom. To enhance the visibility of the 



density perturbation, we subtract a background profile 
n\{z) obtained by averaging the profiles over all mea- 
sured delay times. Our signal 6rii(z, t) = ri\{z, t) — fii(z) 
is finally normalized to the maximum observed density 

^1, max- 

The key point for the detection of second sound is the 
coupling [22, 23] between temperature and density vari- 
ations, which occurs in a systems exhibiting thermal ex- 
pansion. The relevant isobaric thermal expansion coef- 
ficient can be obtained from the EOS and, for our ex- 
perimental conditions, is found to be sufficiently large to 
facilitate the observation of a local temperature increase 
as a dip in the density profiles (Methods). 

First sound clearly shows up in Fig. lb. The initially 
induced hump splits into two density peaks (bright), 
which symmetrically propagate outward at an almost 
constant speed, penetrate into the region where there is 
no superfluid (see dashed lines for the superfluid-normal 
boundary) , and finally fade out in the outer region of the 
cloud. For longer times, we observe (not shown) a weak 
collective breathing oscillation to be excited. 

The excitation with our local heating scheme leads to a 
strikingly different picture, as shown in Fig. lc. The two 
density dips (dark) propagate much slower than the first- 
sound signals. They further slow down when approaching 
the superfluid boundary (dashed lines) and finally dis- 
appear without penetrating into non-superfluid region. 
This behavior is our 'smoking gun' of second sound. 

To extract the two sound speeds from the differential 
profiles 5ni(z,t), we determine the positions of the den- 
sity dips or peaks with Gaussian fit functions. The cor- 
responding time dependence is shown in Fig. 2a, as ex- 
tracted from the profiles in Fig. 1. We now use a third- 
order polynomial to globally fit the time-dependent po- 
sitions (solid lines in Fig. 2a). The sound speeds are 
then obtained as time derivatives of the fit curves and 
displayed as solid lines in Fig. 2b. As this procedure by 
design produces a smooth curve and does not provide suf- 
ficient insight into the uncertainties, we also apply a dif- 
ferent procedure to analyze the same data. We consider 
smaller subsets of adjacent points and extract the local 
speeds by second-order polynomial fits. Corresponding 
results are shown in Fig. 2 by the filled symbols. 

The fact that the axial harmonic confinement intro- 
duces a z-dependence of the linear density m allows us 
to determine the temperature dependence of the sound 
speeds without changing the global temperature T of the 
trapped sample. The key is to define a z-dependent Fermi 
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(Methods) as the natural local 



temperature scale. The corresponding reduced tempera- 
ture T/Tp D has its minimum at the trap center (z = 0) 
and increases with z. The superfluid phase-transition is 
crossed when the critical temperature T c = 0.214(16)T^ D 
is reached. 

In Fig. 3, we show the temperature dependence of the 
two speeds of sound, normalized to the local Fermi speed 
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FIG. 2. Extracting the sound speeds, a, The positions 
of the propagating pulses are shown as a function of time. 
The data points (red and blue symbols for first and second 
sound, respectively) result from individual fits to the pulses 
observed at fixed delay times, and the solid lines represent 
third-order polynomial fits to the time-dependent behavior, 
b, The sound speeds are obtained as derivatives of the fit 
curves (solid lines) and alternatively by analyzing subsets of 
nine adjacent profiles (data points). 



v]P = y2kBT^P/m, where m is the atomic mass and 
ks is Boltzmann's constant. The symbols correspond to 
the data displayed in Fig. 2. The solid lines are derived 
in the same way from the corresponding fit curves. To 
get additional information on the confidence level of our 
results, we have analyzed a number of data sets taken 
under similar conditions as the ones in Fig. 1. The regions 
shaded in pink and light-blue display the maximum range 
of variations considering all our different data sets. 

Our interpretation of the experimental results relies 
on an effective one-dimensional approach that allows to 
solve Landau's two-fluid hydrodynamic equations for a 
highly elongated system [25, 26]. The basic assumptions 
arc a thermal equilibrium in the radial direction and suf- 
ficient shear viscosity to establish a flow field that is in- 
dependent of the radial position. Within this theoretical 
framework and under the local density approximation, 
effective ID thermodynamic quantities can be defined by 
integration over the transverse degrees of freedom, such 
that a thermodynamic quantity q yields a ID counter- 
part qi = 2ir J Q qr dr. We can express the normalized 
speeds of first and second sound as 

I 7 Pl (1) 

and 

"2 / T sj n~T . . 

v]P V 2k B T}Pc p in nl > [Z) 

where P\ denotes the ID pressure (unit of a force), s\ = 
Sx/ni is the entropy per particle, and c p \ = T(dsi/dT) p i 
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FIG. 3. Normalized sound speeds and the ID super- 
fluid fraction, a, Speeds of first and second sound, nor- 
malized to the local Fermi speed and plotted as a function 
of the reduced temperature. The data points and the solid 
lines refer to the data set of Fig. 1, following different meth- 
ods to analyze the raw data (see text). The shaded regions 
indicate the maximum range of variations from analyzing dif- 
ferent data sets. The dashed curve is a prediction based on 
Eq. 1 and the EOS from ref. [17]. The dotted horizontal line 
is the corresponding zero-temperature limit [24] for the speed 
of first sound, b, Temperature dependence of the ID super- 
fluid fraction n a i/ni, with symbols, solid line, and shaded 
uncertainty range corresponding to panel a. In both panels, 
the grey shaded area indicates the uncertainty range of the 
superfluid phase transition according to ref. [17]. 



is the isobaric heat capacity per particle. These ther- 
modynamic quantities can be calculated from the EOS 
as functions of the reduced temperature T/Tp D , as we 
discuss in detail in ref. [26]. The quantities n s i and 
n n x = tii —n-si, which represent the linear number densi- 
ties of the superfluid and the normal component, cannot 
be determined from the known EOS. 

The speed of first sound provides us with an impor- 
tant benchmark for our experimental method and the 
interpretation of the measurements in the ID theoretical 
framework. The experimental results (red symbols and 
upper solid line in Fig. 3a) are in excellent agreement 
with the calculation (dashed line) based on Eq. 1 and 
the EOS from ref. [17], which is a further confirmation 
for the validity of our theoretical approach in addition to 
the recent measurements of the T-dependent frequencies 
of higher-nodal collective modes [18]. 

The measured speed of second sound (blue line and 
data points in Fig. 3a) is observed to decrease with tern- 
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FIG. 4. Superfluid fraction for the homogeneous 
case. The data points and the corresponding uncertainty 
range (shaded region) show the superfluid fraction for a uni- 
form resonantly interacting Fermi gas versus T/T c as recon- 
structed from its ID counterpart in Fig. 3b (Methods). The 
two horizontal error bars indicate the systematic uncertainties 
resulting from the limited knowledge of the critical tempera- 
ture T c . For comparison, we show the fraction for helium II 
(solid line) as measured in ref. [27] and the textbook expres- 
sion 1 — (T/T c ) 3//2 for the Bose-Einstein condensed fraction 
of the ideal Bose gas (dashed line). 



perature, in contrast to first sound. The general behavior 
fits very well to the qualitative predictions in ref. [25] . We 
can now extract the ID superfluid fraction n s i/ni , which 
is the unknown quantity in Eq. 2. The result is presented 
in Fig. 3b, where n a \jn\ shows a smooth increase with 
decreasing T below the critical temperature. 

Wc finally reconstruct (Methods) the temperature de- 
pendence of the superfluid fraction n s /n for the homoge- 
neous 3D case, which has been an inaccessible quantity 
so far. The results, shown in Fig. 4, turn out to be rather 
close to the case of liquid helium II (solid line) [27]. In 
particular, the superfluid part is practically 100% below 
0.6 T c . This behavior is quite different from the one ex- 
hibited by a weakly interacting Bose gas, whose super- 
fluid fraction lies significantly below the data of Fig. 4 
and is well approximated by the condensate fraction of 
the ideal Bose gas (dashed line). In strongly interacting 
quantum fluids, the superfluid and the condensate frac- 
tions simultaneously appear at the phase transition, but 
they exhibit quite different temperature dependencies be- 
low T c . Our experimental results provide a new bench- 
mark for advancing theoretical approaches to calculate 
the superfluid fraction, which is a challenging problem in 
quantum many-body physics. 

From the fundamental point of view, the experimen- 
tally determined superfluid fraction represents a so far 
missing thermodynamic function, which contains infor- 
mation on the spectrum of elementary excitations and 
completes the description of the superfluid in terms of 
universal thermodynamics. With respect to future ap- 



plications, the creation of second sound represents a 
paradigm for controlling the relative motion of a su- 
perfluid with respect to the normal component. This 
may find applications in other situations and geome- 
tries, where resonantly interacting Fermi gases are used 
as model systems for exploring dynamical and transport 
phenomena [28-30]. 

We thank E. R. Sanchez Guajardo for his contribu- 
tions in the early stage of this work. The Innsbruck 
team acknowledges support from the Austrian Science 
Fund (FWF) within SFB FoQuS (project No. F4004- 
N16). The Trento team acknowledges support from the 
European Research Council through the project QGBE. 

METHODS 

Trapping potential. In our hybrid trap [31], 
the tight radial confinement with a trapping frequency 
av/(27r) = 539(2) Hz is provided by an infrared laser 
beam (wavelength 1075 nm, power 120 mW, waist 39/im). 
The much weaker axial confinement with a trapping fre- 
quency oj z /(2ir) = 22.46(7) Hz results from the curva- 
ture of the applied magnetic field. The frequency ra- 
tio u> r /uj z w 24 corresponds to the aspect ratio of the 
trapped cloud. 

Sound excitation and detection. The green laser 
beam (wavelength 532 nm, power < 20 mW) used for ex- 
citation is focused to a waist in the range of 25-35 fim. 
To excite first sound, we suddenly turn it on to introduce 
a repulsive potential hill with a height of about 10% of 
the cloud's chemical potential in the trap center. For the 
case of second sound, the repulsive beam is permanently 
on during the preparation of the quantum gas with a bar- 
rier height of about 15% of the cloud's chemical potential. 
The excitation is then induced by a burst of 8 cycles of 
sinusoidal power modulation. The modulation frequency 
is set to 1720 Hz, but the scheme is rather robust and was 
found to work best in a frequency range between roughly 
2 and 3.5 times the radial trapping frequency. The burst 
envelope is rectangular and the modulated barrier has 
a pcak-to-pcak amplitude of 30% of the cloud's chemi- 
cal potential. After the burst, the power is set back to 
its initial constant value. On the time scale of the axial 
motion, the time-averaged power of the green beam is 
constant, which avoids direct excitation of first sound. 

For the detection of the propagating second-sound sig- 
nal, the corresponding density dip is essential. The ID 
formulation of the EOS in ref. [26] allows us to relate 
the observed depth to the relative temperature change. 
In the temperature range of our experiments, the cor- 
responding thermodynamic coefficient (8n\ / n{) / (5T / T) 
takes values [26] between —0.4 and —0.6, which means 
that the typical 3% relative depth of the density dip 
roughly corresponds to a local temperature increase of 
about 6%. 
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ID Fermi temperature and critical tempera- 
ture. Three different definitions for Fermi tempera- 
tures are related to natural energy scales of our trap- 
ping geometry. The homogeneous case with a 3D num- 
ber density n (including both spin states) is given by 
k B T F = (37r 2 ) 2 / 3 ^ri 2 / 3 . Within the local density ap- 
proximation the corresponding Fermi energy for N atoms 
in a three-dimensional harmonic potential is given by 
kBTp ap — h(3N td'^ujz) 1 ^ 3 , commonly used to describe 
the global situation of a 3D trap. The Fermi energy of 
cylindrically confined cloud in the center of a 2D har- 
monic trap, which we refer to as 'ID Fermi temperature', 
follows from k B T F D = (l|L)a/B(fi Wr )4/6(^)i/S j wne re 
the relevant density is the linear number density ri\. For 
the additional axial confinement in our trap geometry, n\ 
and thus T^P become z-dependent quantities. 

For specifying the critical temperature T c in units of 
the relevant Fermi temperature, we also distinguish be- 
tween the three different situations of a homogeneous 
system, a 3D harmonic trap, and a cylinder with 2D 
harmonic confinement. For the homogeneous system, 
T c = 0.167(13) Tp was measured in ref. [17]. Based on 
the local density approximation and the experimentally 
determined EOS, this result can be translated into corre- 
sponding conditions for the other two situations. While 
for the 3D trap, T c = 0.223(15)x£ rap is relevant for the 
occurrence of the phase transition at the center of the 
trap, the condition T c = 0.214(16)7^° applies to the 2D 
confined case. In our highly elongated trap geometry 
with weak axial confinement, Tp^ and thus T c become 
z-dependent. For a fixed global temperature T, the con- 
dition T < T c (z) then determines the axial range, where 
a supcrfluid exists (see illustration in Fig. la and dashed 
lines in lb and lc). 

Reconstruction of superfiuid fraction. Within the 
framework of universal thermodynamics [13], the number 
density of a uniform, resonantly interacting Fermi gas can 
be expressed in terms of a dimcnsionlcss universal func- 
tion f n (x) by n(x,T) = X T 3 f n (x). Here A T = {^^f 2 



is the thermal dc Broglie wavelength, and the dimen- 
sionless parameter x = fi/k B T gives the ratio between 
the chemical potential y, and the thermal energy k B T, 
with a unique correspondence existing between x and 
T/Tf- The function f n {x) is known from measurements 
of the EOS [14-17]. The supcrfluid density can be ex- 
pressed in an analogous way as n s = X T 3 f ns (x), intro- 
ducing a corresponding universal function f ns (x), which 
is to be extracted from our measurements. Using the 
local density approximation, one can show for a system 
with radial harmonic confinement that [16] tii(xq,T) = 
-^ZP^yjf- Jl^ fn(x)dx, where xq represents the value of 
x on the trap axis. Analogously, the ID supcrfluid 



only depends on xq. From our experimentally deter- 
mined n s i/ni, we readily obtain the supcrfluid fraction 
of a uniform gas using the relation — — f TlB ( x °' ) 

^r_if n ( X )dx 



fn(x ) 



1 d 

fn(xo) dx 



density is given by n sl (x ,T) = f„ s (x)dx. 
One thus easily sees that the ID superfiuid fraction 
is given by n s \/m = f ns (x)dx/ f* ^ f n (x)dx and 
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